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THE PROFILE OF BOUNDARY GRADIENT BLOW-UP 
FOR THE DIFFUSIVE HAMILTON-JACOBI EQUATION 

ALESSIO PORRETTA AND PHILIPPE SOUPLET 


Abstract. We consider the diffusive Hamilton-Jacobi equation 

ut — Au = |Vm|’’, 

with Dirichlet boundary conditions in two space dimensions, which is a typical model-case in the 
theory of parabolic PDEs and also arises in the KPZ model of growing interfaces. For p > 2, 
solutions may develop gradient singularities on the boundary in finite time, and examples of single¬ 
point gradient blowup on the boundary are known, but the space-profile in the tangential direction 
has remained a completely open problem. In the parameter range 2 < p < 3, for the case of a flat 
boundary and an isolated singularity at the origin, we give an answer to this question, obtaining 
the precise final asymptotic profile, under the form 

+ \ as ( 1 , 1 /)-s-(0,0). 

Interestingly, this result displays a new phenomenon of strong anisotropy of the profile, quite dif¬ 
ferent from what is observed in other blowup problems for nonlinear parabolic equations, with the 
exponents l/(p — 1) in the normal direction y and 2/(p — 2) in the tangential direction x. Further¬ 
more, the tangential profile violates the (self-similar) scale invariance of the equation, whereas the 
normal profile remains self-similar. 


Contents 

1. Introduction [1] 

1.1. B ackground [1] 

1.2. Main result: final gradient blowup profile near an isolated boundary singularity [2] 

1.3. Discussion and remarks 0] 

1.4. Existence of single-point gradient blow-up solutions [5] 

1.5. Ideas of proof [6] 

2. Preliminary properties [7] 

3. Proof of main result: the lower estimate m 

4. Regularizing barriers [16] 

5. Proof of main result: the upper estimate [T8| 

6. A heuristic explanation of the singularity exponents through quasi-stationary 

approximation [25| 

References [26| 


1. Introduction 

1.1. Background. This article is devoted to the qualitative study of solutions of the diffusive 
Hamilton-Jacobi equation 

Ut-/S.u=\Vu\P. (1.1) 

2010 Mathematics Subject Classification. 35K55, 35B40, 35B44, 82C24. 

Key words and phrases. Diffusive Hamilton-Jacobi equations, KPZ model, gradient blow-up, isolated boundary 
singularities, blow-up profile, tangential profile, anisotropic singularity. 

1 



2 


PORRETTA AND SOUPLET 


Beside being the viscosity approximation of Hamilton-Jacobi type equations from stochastic control 
theory [28], equation (1.1) is involved in certain physical models, for instance of ballistic deposition 
processes, were it describes the evolution of the profile of a growing interface. It is actually the 
deterministic version of the well-known Kardar-Parisi-Zhang (KPZ) equation (see |23j and cf. also 
Krug and Spohn |25| ). In its stochastic version, it has undergone spectacular development recently 
with the work of M. Hairer |19) . Finally, equation (1.1) is a typical model-case in the theory 


of parabolic PDFs. Indeed it is the simplest example of a parabolic equation with a nonlinearity 
depending on the first-order spatial derivatives of u. As such, it is important to study its qualitative 
properties. 


Equation (1.1) has been intensively studied in the past twenty years, and it is well known that 
two fundamentally different situations occur. If the equation is considered in the whole space 
(with, say, bounded initial conditions), then all solutions exist globally in the classical 
sense and remain bounded in see e.g. lailllHlEelClIITlIIlllMlESI, where the large time 
behavior is investigated in detail). At the opposite, if the equation is posed on a domain with 
homogeneous Dirichlet boundary conditions, then for p > 2 and suitably large initial data, the 
local classical solution develops singularities in finite time. These singularities are of gradient 
blowup type (GBU), the function u itself remaining bounded, and are located on some part of the 
boundary; see e.g. 12llI31[IIllIlE3ia[22lElEaiMlE21IISlEZj and the references therein. 

For the classical blowup problem associated with the nonlinear heat equation 


ut — Am = 


( 1 . 2 ) 


a considerably developed theory is available for the description of the asymptotic profile of the 
solution near a finite time singularity (see |32j and the references therein). In comparison, very 


little is known for equation (1.1). In particular, in the case of an isolated boundary singularity, the 

final blowup profile of Vm in the tangential direction is completely unknown]^ 

1.2. Main result: final gradient blowup profile near an isolated boundary singularity. 

The goal of this paper is to fill this gap by giving a substantial contribution to this question. In 
the range of exponents 2 < p < 3, we will give a sharp description of the final blowup profile of 
Vm near an isolated boundary singularity (in both normal and tangential directions). Since the 
question is quite involved, we shall restrict ourselves to a rather simple setting, but which captures 
the essence of the problem. Namely, we consider the two-dimensional case, where the domain is 
assumed to coincide locally with a half-plane near the point of singularity. To this end, for given 
p > 0, we set 

ujp = {-p, p) X (0, p) C a;+ = Wp n {x > 0}. 

Next we fix some L,T > 0 and put 

uj = u;l, J =ojlI2-, (1-3) 

gT:=a;x(0,r), Ft = (-L, L) x {0} x (0, T). (1.4) 

Definition 1.1. Let L,r > 0 and let u G x (0,T)) be a nonnegative classical solution 

of ( fj.ip in Qt, with u = 0 on Tt- We say that u has an isolated gradient blowup point 
at (0,0, T) if 

limsup I Vm(x, y, t)| = oo (lA) 

and 

Vm is bounded on K X {0,T) for any K CCu}\{{0,0)}. (1.6) 


^It is only known that |Vu(X, t)| < C[dist(X, (see 1^ 15^ I27 |l. which gives an upper estimate in 

the normal direction but provides no information on how the profile is damped away from the point of singularity 
along the boundary. 
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If u has an isolated gradient blowup point at (0,0, T) then we may define the final blowup 
profile of Vu, given by 

Vu(a:, y, T) := lim Vu{x, y, t), for all (x, y) GlJ' \ {(0,0)}. 

t—^T 


Indeed the limit above exists and is finite due to (1.6), as a consequence of standard parabolic 
estimates. Our main result is the following. 


Theorem 1.1. Assume 


2 < p < 3. 


Let L,T > 0, let u G x (0, T)) be a nonnegative classical solution of 

on Tt- Assume that u has an isolated gradient blowup point at (0,0, T) 
monotonicity condition 


( 1 . 1 ) in Qt, with u = 0 
and that u satisfies the 


xux <0 in Qt- 


(1.7) 


Then there exist constants C'i,C' 2 ,C '3 > 0, p G (0,L) (possibly depending on u) such that, for all 
{x,y) G ([—p, p] X [0,p]) \ {(0,0)}, the final blowup profile satisfies 


-0 


y + C7i|x|2(P-i)/(P-2) -C3<Uyix,y,T)<dp y + C2 |x|2(p-i)/(p-2) 


-0 


+ C3 


( 1 . 8 ) 


where 

j3 = l/{p—l) and dp = (3^. 

In particular, the final profile of the normal derivative on the boundary satisfies 

C4|xr2/(P-2) < Uy{x,0,T) < 

for all 0 < |x| < p and some C 4 , C 5 > 0. Also, for some C > 0 , we have 

u < C, \ux\ < C, for all {x, y) G oj'. 



Fig. 1: The shape of the final profile of u near the origin. 
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1.3. Discussion and remarks. 

(a) Interestingly, this result shows that the GBU profile is strongly anisotropic, i.e. the expo¬ 
nents of the singularity profile in the normal and in the tangential directions are different, respec¬ 
tively l/(p — 1) and 2/(p — 2). Moreover, whereas the exponent of the normal profile obeys the 
natural scaling of the equation, the latter is violated by the tangential profile. Indeed, recall that 
equation 0 is invariant under the group of transformations 

u ux{x,y,t) := X^u{Xx, Xy, X^t) with m = (2 — p)/(p — 1), for all A > 0, 

whose gradient is given by Vux = X^^^p~^^Vu{Xx, Xy, A^t). 

(b) As far as we know, no similar example of anisotropic, isolated blowup singularity is known in 


parabolic problems. For the nonlinear heat equation (1.2), the stable blowup profile at an isolated 
blowup point is known to be isotropic]^ (see [371130l EH US] and the references therein), with 

u{X,T) ~ c(p)|A:|-2/(P-i)|log|A:||-i/(P-i) as A ^ 0. 

Here X G M” with n > 2 and 1 < p < {n + 2)/{n — 2), and this profile occurs for instance for any 
symmetric, radially decreasing solution. 

The case of the linear heat equation with nonlinear boundary conditions 


ut — Art = 0 in H X (0, T), 

— =uP on dn X (0, T) 
av 


(1.9) 


was studied in uni EH [201 EH . Like for ( jl.lj ), this problem involves boundary singularities (however 
u itself blows up). It was recently found in [20l EH that for Q = = {{x,y); y > 0} under 

assumption (1.7), the singularity profile satisfies 

-i/(p-i) 


u{x,y,T) 


y 


X l/(p ^)|logx| 


for y —)■ 0 with |a:| = 0 {y) 
for X —)■ 0 and y = 0. 


A similar result holds in dimension n > 3 if 1 < p < n/{n — 2). Note that this profile is only weakly 


anisotropic (by a logarithmic correction) in comparison with (1.8). 


On the other hand we also observe that, unlike in problems (1.2) and (1.9), the profile that we 


find for (1.1) is given by pure powers, without (e.g. logarithmic) corrections. This situation seems 


typical of type II blow-up problems (see [29] and cf. Remark (c)). 

(c) The exponent 2j{p — 2) appears to be new in this problem. However, it is worth noting that, 
in some cases, the time rate of GBU involves a related exponent. Namely, for monotone in time 
solutions in 1 space dimension, we have [18j : 

||Vn(-,t)||oo~(r-t)-'/(^-"). (1.10) 

However, the question of the time GBU rate is still open in 2 dimensions. Note that the rate ( l.lOj ) 
corresponds to a type H blow-up, in the sense that this rate is more singular than what the natural 
scaling of the equation would suggest (see [Ml ES] for details). A possible heuristic explanation 
of the appearance of the number 2/(p — 2) in this problem, based on ideas of quasi-stationary 
approximation, is given in Section 6. 

(d) It remains an open problem what is the actual tangential singularity exponent for p > 3 

- see Remark 6.1 for details. Actually the lower estimate in remains true for any p > 2 

(cf. Theorem 3.1). As for the upper estimate, for p > 3, our method would allow to obtain an 


estimate of the form in (1.8), with some power, greater than 2 {p — l)/(y — 2), which could be 


explicitly computed in terms of p. However, due to the gap between the upper and lower estimates 


^this concerns blowup at an interior point - actually only the whole space case is considered in these works; 
however no blowup can occur at a boundary point for equation (1.21, at least in a convex domain 
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in this case, we are unable to determine the exponent of the actual prohle. Therefore, and in order 
not to further increase the technicality of the article, we have refrained from expanding on this. 
On the other hand, it might be possible to extend our results to more general (nonflat) domains 
and to higher dimensions, at the expense of further complication. But since the main goal of this 
work is to present a new phenomenon, we have decided to leave this aside. 


(e) Actually, the upper estimate in (1.8) is satisfied by Uy{x,y,t) for all t < T (this is a conse¬ 


quence of the proof, cf. in particular formula (5.32)). 


(f) By direct integration of (1.8) between 0 and y, one easily obtains the corresponding estimate 


for the profile of the function u itself (whose shape is depicted in Fig. 1). In the course of the 
proof of Theorem [13 we also establish additional estimates, of possible independent interest. In 
particular, we show that for any p > 2, there holds 

\ut\ < C, \Ux\ < C\x\, Uxx > —C 
for (x, y, t) close to (0,0, T). Moreover, for 2 < p < 3, we show that 

-C < Uxx{0, y,T) < -c < 0 


for y > 0 small (see Remark 5.1). In particular, since Uxxix,0,T) = 0 for x 7 ^ 0, we see that 
Uxx{--,T) is discontinuous near the origin. 


1.4. Existence of single-point gradient blow-up solutions. In order to obtain solutions sat¬ 
isfying all the assumptions in Theorem o we now recall a result from m concerning the initial¬ 
boundary value problem 


ut — Au = |Vu|^, 
u{x, t) = 0 , 
u{x,0) = uo{x), 

Here, it is assumed that 

H C is a C^+^-smooth bounded domain, uq G C^(H) with uq > 0 and Jto|an = 0. 


X G H, t > 0, 

X G dQ, t > 0, 
X G H. 


( 1 . 11 ) 
( 1 . 12 ) 

(1.13) 

(1.14) 

It follows from m Theorem 10, p. 206] that problem ( |1.11[ )-( [1.13[ ) admits a unique maximal, 
nonnegative classical solution u G x (0,T)) n C'^''^(H x [0,T)), where T = T{uo) is the 

maximal existence time. Also, by the maximum principle, we immediately have 

||jj(i)||oo < ||wo||oo, 0 < t < T. 

On the other hand, by the Bernstein-type estimate in |36j . we know that 

\Vu{t,X)\ < C[dist(A,dH)]-i/(P-i), in H x (0,r), 
so that in particular GBU can take place only on cAI. The following result was proved in m- 


Theorem A. Assume and 

Q, and uq are symmetric with respect to the line x = 0, (1-15) 

H coincides locally near 0 with the half-plane {y > 0 } and is convex in the x-direction, (1.16) 
xuq^x < 0 in Q. (1-17) 

If Uq is suitably concentrated near the origin (see Remark M below), then the solution of {1.11^- 


(1.13) satisfies 


T = T{uq) < 00 and Vu blows up only at the origin 


(i.e. (1.5) is true and Vu is bounded on K x (0, T) for any K CC H \ {(0,0 )}). 


Also we note that, as a consequence of the assumptions of Theorem A, we have xux < 0 in 

n X (o,r). 
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Remark 1.1. As an example of data “suitably concentrated near the origin” for Theorem A, the 
following was given in m- 


uo{x) = Ce’"ip(^£ + (y - (1-18) 

where k = {p — 2)/{p — 1), C > Co{p) > 0, e > 0 is sufficiently small, and p G C'°°([0,(X))) satisfies 
p' <0, ¥5(s) = 1, s < 1/4, (^(s) = 0, s > 1/2. (1.19) 


Note that these functions have support concentrated near the boundary point (0, 0) (and large 
derivatives for e small). 


1.5. Ideas of proof. The proof of Theorem 1.1 is long and technical, and it requires to combine 
many ingredients. Let us briefly describe the main ideas. 


To establish the lower estimate, we start with an estimation of the normal derivative on the 
boundary, which is obtained in three steps (see Fig. 2 in Section 3): we start from the vertical 
line {x = 0}, where the precise final profile follows rather easily from ODE arguments. We then 
extend the lower estimate to the region above the curve y = by using a lower bound of 

Uxx along horizontal segments. The extension to the region below the curve y = is then 

achieved by means of a boundary Harnack-type estimate in suitable boxes connecting this curve to 
the boundary {y = 0}. Once Uy is estimated from below on the boundary, the full lower estimate 
of Uy is obtained by suitable integration along vertical lines, plus some horizontal averaging made 
possible by an estimate of the mixed derivative Uxy 

As for the proof of the upper estimate, it combines two ingredients. The first one is an auxiliary 
function of the form 


J{x,y,t) = Ux + kx{l + y)y 


with suitable parameters k,q > 0, which is shown to be nonpositive via the maximum principle. 
The integration of the inequality J < 0 along horizontal lines then yields a sharp upper estimate of 
Holder type for u. The second ingredient is a family of suitable regularizing barriers, which allow 
us to improve the Holder estimate of rt to a pointwise upper bound of Uy on the boundary. We note 
that rougher versions of both ingredients were used in 1211, in order to prove single-point GBUj^ 
However, these ideas need to be considerably refined in order to obtain the sharp tangential GBU 
profile. In particular, the derivation of the parabolic inequality satisfied by J requires a delicate 
analysis in terms of the auxiliary quantities 


i = and Q = y{uyf ^ 

This latter step turns out to require the restriction p < 3 and leaves open the question whether the 
upper estimate remains true for p > 3 as well (see Remark 

The rest of the paper is devoted to the proof of Theorem o Some preliminary properties, 
mostly based on the maximum principle, are given in Section 2. The lower estimate is established 
in Section 3. In Section 4, we construct the regularizing barriers. In Section 5, the analysis of the 
parabolic inequality for the function J is carried out, and the proof of the upper estimate is then 
completed. Finally, a possible heuristic explanation of the appearance of the number 2/(p — 2) in 
this problem, based on ideas of quasi-stationary approximation, is given in Section 6. 


^The function J in m was itself motivated by a device from m, where a func tion of the form J{r,t) = 
+ er"‘u'^ was introduced to study the blowup of radial solutions of equation (1.2 I. 
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2. Preliminary properties 


In the following propositions, we state a number of useful bounds and properties of the solution, 
which will be used in the proof of the main result Theorem 1.1. All the proofs will be given after 
the statements. Here and in the rest of the paper, letters such as (7, Ci, c,... will denote possibly 
different positive constants, whose dependence will be specified only when necessary. 

We start with some simple bounds, which follow rather easily from the maximum principle. Let 
us recall that u^oj',Q t and Pt are defined in (1.3)-(1.4). 

Proposition 2.1. Assume p > 2, let L,T > 0 and let u G x (0,r)) be a nonnegative 

classical solution of (1.1) in Qt, with u = 0 on Pt^. Assume that u has an isolated gradient blowup 
point at (0, 0, T). 

(i) Then u extends to a function 

ueC^’^iQ), withQ := {oJ'x[T/2,T])\ {{0,0, T)}. (2.1) 

(This extension will still be denoted by u without risk of confusion.) 

(a) There exists a constant C > 0 (possibly depending on the solution u), such that u satisfies 
the following bounds in Q: 


u 


\ut\ < c 

(2.2) 

Uy > —C 

(2.3) 

Uxx > —C. 

(2.4) 

Y1 ^ 1 ^ 1 

(2.5) 


If, moreover, u satisfies (l.I), then we have 
in Q. 

We next show that the gradient blowup does occur in a pointwise sense: Uy becomes uniformly 
large near the blow-up time and the origin. 

Proposition 2.2. Assume p > 2, let L,T > 0 and let u G x (0,r)) be a nonnegative 

classical solution of (1.1) in Qt, with u = 0 on P^. Assume that u has an isolated gradient blowup 
point at (0, 0, T) and that u satisfies the monotonicity condition ( 73 . Then we have 

lim Uy{x,y,t) =+CO. 


t^T 
(x,y)^{0,0) 


( 2 . 6 ) 


As a consequence of Proposition |2.2[ there exists 0 < po < min(L/2,r/2) such that 

uy>i in (ojpo X [r -/ 9 o,r])\{(o,o,r)}. (2.7) 

We now give upper bounds which essentially follow by integrating in the vertical direction. 

Proposition 2.3. Under the assumptions of Proposition there exist a constant C > 0 and 
p G (0,/Oo) (possibly depending on u), such that the solution u satisfies 

Uy{x,y,t) < [{uy)^~P{x,0,t) + {p- 1 ) 2 /]"^ + 2Cy in ujp x [T - p,T). 


In particular, we have 


where 


Uy{x,y,t) <dpy ^ + 2Cy 
u{x,y,t) < Cpy^~^ + Cy‘^ 


in ujp x\T — p, T) 
in ujp x\T — p, T), 


(3 = 


1 ’ 

P — 1 


dp = 133 , Cp = {I - (3) ^dp. 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 
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Our next result shows that similar lower bounds are true at x = 0 (of course they cannot be true 
for X 7 ^ 0 in view of the profile eventually found in (1.8)). 


Proposition 2.4. Under the assumptions of Proposition \2^ there exist a constant Oi > 0 and 
p G (0,/Oo) (possibly depending on u) such that 


/{0,y,t) > [{uy)^ P{0,0,t) + ip-'i-)y] ^-Cl. 0<y<p, T-p<t<T. 


Moreover, we have 
and 


uy{0, y,T)>dpy ^ - Cl, Q <y < p 


( 2 . 12 ) 


(2.13) 


tt(0, y, T) > Cpy^ ^ - Ciy, 0 < y < p. (2.14) 

The following relationship between second order derivates, whose proof is rather delicate, will 
play an important role to establish the lower pointwise estimates in (1.8). 

Proposition 2.5. Under the assumptions of Proposition\2^ for any p > 0, there exists a constant 
Cp > 0 (possibly depending on u), such that the solution u satisfies 


Uxy<puxx + Cp in X [T/2,T]. 
We now turn to the proofs of the results that we have just stated. 


(2.15) 


Proof of Proposition 2.1 


• Proof of (2.2) and (2.3). Set z = ut or z = Uy. 
it satisfies 

zt — Az = A 

where A = p\Vu\^~‘^Vu. 


Property (2.1) is a consequence of standard parabolic estimates. 

Then z G C^’‘^{Q) by parabolic regularity and 


Vz, 


(2.16) 


Since u* = 0 on T^, using (2.1), we see that the supremum of \ut\ on the parabolic boundary of 
u' X [T/2,T) is finite. Denoting this supremum by C, the maximum principle then guarantees 

\ut\ < C in cj'X [r/2, r), 


which implies (2.2) in view of (2.1). 

We can apply a similar reasoning to Uy. Since u> 0 and u = 0 on T^, we have tty > 0 on T^. By 


(2.1), we see that the infimum of Uy on the parabolic boundary of uj' x [T /2, T) is finite. Denoting 
this infimum by — C, the maximum principle then guarantees 


Uy>—C in cj' X [r/2, T), 

which implies 
• Proof of 

Zt-= p[\Vu\P-‘^Vu-Vux\^ 

= A-VZ + p\Vuf’-‘^\Vux\'^ + p{p - 2)\Vu(P-^{Vu ■ Vuxf (2.17) 

>A-VZ. 


(2.3) in view of (2.1). 


(2.4). The function Z := Uxx G by parabolic regularity and it satisfies 


Since Z = 0 on Tt^, using (2.1) we see that the infimum of Z on the parabolic boundary of 
uj' X [T/2,T) is finite. It then follows from the maximum principle that Uxx P —C in uj' x [T/2,T), 
which implies (2.4) in view of (2.1). 


Proof of (2.5). As a consequence of (|1.7[), we have 


Ux{0,yU) = 0 for all {y,t) G ([0,L/2] x [T/2,r]) \ {(0,T)}. 


(2.18) 
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Consequently, we get 


Ux{x,y,t)= / Uxx{t,s,y)ds> -Cx in (w'x [r/2, T)) n {x > 0}, 
Jo 


and a similar estimate for x < 0. This implies (2.5). 

In view of the proofs of Propositions |2.2^2.41 we prepare the following Lemma. 
Lemma 2.6. Under the assumptions of Proposition we have 

limsup%(0,0,t) = +00 


□ 


t-s-T 


and 


Uy{0, y, T) > dpy ^ - Cl, 0 <y < 


Proof. As a consequence of (1.7) and u = 0 on P-r, we have 


hence 


xuyx{x, 0, t) < 0 for 0 < |x| < L/2 and f G [T/2, T], 


Uy{0,0,t) = sup Uy{x,0,t) for t G [T/2,T). 
|x|<L/2 


(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 


On the other hand, by (2.16) and the maximum principle, for each r G (T/2,T), the maximum of 
Uy on Q(. := oj' x (T/2,r) is attained on the parabolic bonndary dpQf of Qf. Moreover, by (1.6), 
we have 


Mo = sup Uy < oo, where P' := dpQ'rp \ ([—L/2, L/2] x {0} x [r/2,T)). 


Therefore, 


supuj; < max( Mq, sup Uy(0,0,t)). 

Q'r ^ ie[T/2,r] J 


( 2 . 22 ) 


By our assumption (1.5), the LHS goes to oo as r —)• T and (2.19) follows. 

Let us now prove ( 2.20| ). By (1.7), we have 

Uxx{0,y,t) < 0 for all {y,t) G (0,L/2] x (r/2,r]. 

Also, we know from Proposition 2.1 that Uy > —C hence \uy\ < Uy + 2C and that \ut\ < C in 
u' X (r/2,r). Set Cl := 2C + C^. For {y,t) G (0,L/2) x {T/2,T), using (2.18) and (2.22) it 
follows that 

-Uyy{0,y,t) = [uxx + |Vu|P - ut\{0,y,t) < \uy{0,y,t)\P + C < {uy{0,y,t) + CiY- 
Observing that Uy{0,y,t) + Ci > 0 and integrating in y, we obtain 

(uj/(0, y, t) + Cif-~^ < (uy(0,0, t) + Cif~P + {p- l)y- (2.23) 

By (2.1) and (2.19), we deduce that 

{uy{0,y,T) + Ci)^~P < liminf(uy(0,0,t) + Ci)^~p + {p - l)y = {p - l)y for y G (0, L/2), 
which yields ( |2.20| ). □ 


Proof of Proposition \2^ Assume by contradiction that there exist a constant K > 0 and a sequence 
{tn,Xn,yn) SUch that 


ixn,ynitn) t (0,0,L) and Uy(^Xji,ymtn) ^ at. 


By (2.2) and (2.4), we have 


^yy — ^yy [Vrij^ — ut Uxx — 2(17, 


(2.24) 
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hence in particular Uyy < 2C. Fix any y G (0, L/2). For n large enough, we have 0 < yn < y, hence 

Uy (^Xn, y^tn) ^ UyiXyi^ y-n: in') ‘2C{y yn) ^ K -|- CL. 

Letting n —>• oo and using we get 

Uy(0, y, T) < iL + CL. 

This is in contradiction with (2.20). □ 


Proof of Proposition 2.3. Fix any x G {—Ll2,Ll2) and t G {T/2,T). By (2.7), there exists 
p G (0,/9o) such that 

h{y) := Uy{x, y, t) - 2Cy > 0, for y G (0, p). 

By (2.24), the function h satisfies 

h' + W = Uyy — 2C + (Uy — 2C y)^ < Uyy — 2C + {Uy)^ < 0. 

By integration, we obtain 

h{y) < [/i^"^(0) + {p- l)y] for y G (0, p), 

hence (2.8) and in particular (|2.9[). Property (|2.10[) follows by further integration. □ 


Proof of Proposition 2.4 Estimate (2.12) is an immediate consequence of (2.23). As for (2.13), it 
was already proved in Lemma 2.6 Finally, (2.14) follows from (2.13) by integration. □ 


We finally prove Proposition 2.5 


Proof of Proposition \2.f^ In view of estimate (2.4), there is no loss of generality if we only consider 
r/ < 1. First we recall from ( |2.17 ) that Uxx satisfies 

{uxx)t — ^XLxx — A • Vuxx > 0 , A = p\Vuf’~‘^Vu . 

We compute the same equation for Uxy, and we get 

(Vm • Vux)iyu ■ Vuy)' 


(^^xy)t ^Uxy A • ^Uxy — ^Uy • ^ Ux T (p 2)”' 


|Vm|2 

= \Vu\^Vuy ■ Vux + {p — 2)(Vu ■ Vux){Vu ■ Vuy) 


(2.25) 


Notice that this is justified close enough to the singularity, due to > 0 (cf. ( 2.27[ ) below) and 
parabolic regularity. Now, given y < 1, we consider the function 

-2^ ■ — '^xy V ^xx 


and the operator C{z) = zt — Az — A ■ Vz. On account of (2.25) we have 

C{z) < F in Qt 

and 

F := \Vu\‘^Vuy ■ Vux + (p — 2){Vu ■ Vux){Vu ■ Vuy) . 

We analyze now the sign of F at large values of 2 ;. First of all, we develop each component of 
the scalar products and we find 

P 

r = Uxy {uxx{{p — ^){Ux) + {Uy) ) F Uyy{{Ux) + {p — l){Uy) )} /r, 

(2.26) 


p\Vu\P' 


T (p 2)UyUxUyyUxX T (P 2)UxUyU, 


xy 


Due to Proposition 2.2, we may choose r > 0 small so that Uy is sufficiently large in Q'j. := 
(0,r)2 X (T — r,T). In particular, we may assume that 

Uy ^ 0, ^yy F 0, nxx F ^yy iP Qr^ (2.27) 
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the last two inequalities coming from the equation Uxx + Uyy = — |Vu|^ + ut together with the 


bounds (2.2) and (2.4). 


In addition, since C wj x (0,T), we have Ua; < 0 and therefore the last term in (2.26) is 
nonpositive. Dropping this term we get 

P 

134 < Uxy {uxx{{P ~ l)('Wx) + (Uy) ) + '^yyii'^x) + (P ~ l)('^y) )} 


p\Vu\P~ 


(P ‘^'}'^yUxUyyUx 


Now, since Uxx > —C (cf. (2.4)), at any point where z > My := pC we have Uxy > My+rjUxx > 0, 
e e 

F 


and since Uxx < ~'^yy ''^6 estimate 


p\Vu\P~^ 


_4 ^ (p 2) UxyUyyi^i^Uy^ {Fx') ) 4“ '^y'^x'^yy'^xx 


— (p (ua:) ) 4“ UyUxUx 


(2.28) 


^y fxx fxxx 

We conclude by noticing that the right hand side of (|2.28|) is negative if Uy is large enough. Indeed, 


by Proposition 2.2 and (2.5), we may chose r so that 


Uy ^ \\UoL 
7] 


in Q'^. 


Then, at any point of Q'^ where z > M„, we have 


'^xy ^ My + pUxx — 4“ C*) > {Uxx 4“ C*) ^ ^ Uxx ; 


U,, 


hence 


Uxy(^{Uy^ {Ux) ) 4“ UyllxUxx ^ (^^xyiUy^ 4“ UyUxUxx — "^yi^c^^xyUy 4“ UxUxx^ ^ 0. 


So from (2.28) we get T < 0 at any point of such that z > r]C. On the other hand, considering 


the parabolic boundary of we have z < r/C at {x = 0} due to (2.18) and (2.4). At {y = 0}, we 


have Uxy — rjUxx = Uxy < 0 by (2.21). On the rest of the lateral boundary, as well as at t = T — r, the 


function is bounded by some constant Cy. By the maximum principle applied to £, we deduce that 
z < max(r/C, Cy) in Q'^. Therefore, the bound (2.15) is proved in In view of the regularity of u 
outside of (0,0, T), the bound can of course be extended to ivj) x [r/2,r] up to an extra uniform 
constant. □ 

3. Proof of main result: the lower estimate 

In this section, we shall prove the following Theorem, which is valid for any p > 2, and in 
particular implies the lower estimate in Theorem o 


Theorem 3.1. Assume p > 2, let L,T > 0 and recall the notation in {1.3)-(1.4)- Let u G 
C'^’^(tJ X (0,T)) he a nonnegative classical solution of (1.1) in Qt, with u = 0 on Ty. Assume that 
u has an isolated gradient blowup point at (0,0, T) and that u satisfies the monotonicity condition 

xux <0 in Qt- (3T) 

Then there exist constants Ci,C 2 > 0, p ^ (Oj-^) (possibly depending on u) such that, for all 
(x, y) GuJp\ {(0, 0)}, the final blowup profile satisfies 


where 


Uy{x,y,T) > dp y + Ci\x\‘^^P 

fi = l/{p—l) and dp = 


1-/3 


-^2, 


(3.2) 
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In particular, the final profile of the normal derivative on the boundary satisfies 

nj,(x,o,r) >C3 |x|-2/(p-2), 
for all 0 < |x| < /> and some C3 > 0. 

In the rest of this section we denote the final prohle at the blow-up time by 

u:=u(-,r)GC2(w'\{(0,0)}). 


(3.3) 


Theorem 3.1 is proved in two steps. First, in Lemma[3.2[ we establish the estimate of the normal 


derivative on the boundary (i.e. (3.3)). To do so, the idea is as follows (see fig. 2 below): we start 
from the vertical line {x = 0}, where the precise lower bound of the final profile v is already known 
thanks to (2.14). We then extend the lower estimate of v to the region above the curve 


So = {(x,y) = (3.4) 

which plays an important role in our arguments. This relies on the lower bound of Uxx in Proposi¬ 
tion used along horizontal segments. This is not sufficient since the region S_|_ does not touch 


the boundary {y = 0}. However, the extension to the region S_ below the curve (3.4) can then 
be achieved by using a Harnack-type estimate in suitable boxes connecting the curve Sq to the 
boundary {y = 0}, in terms of the distance to the boundary. 

Finally, once the normal derivative is estimated on the boundary, the full lower estimate of Uy is 


obtained (cf. Lemma 3.4) by suitable integration along vertical lines, plus some horizontal averaging 
made possible by the estimate of the mixed derivative Uxy given in Proposition 2.5 


Lemma 3.2. Under the assumptions of Theorem 3.1, there exist constants cq > 0 and p G (0,L) 
such that we have 

Vy{x, 0) > for 0 < X < p. (3.5) 


For the proof of Lemma |3.2[ we shall use a well-known quantitative version of the Hopf Lemma 
(or boundary Harnack inequality) [H], which we state in a suitably scale invariant form. 

Lemma 3.3. Let Di be a domain such that 

(-1,1) X (0,2) CZ4i C (-2,2) X (0,2). 

For any (xo,yo) £ <^i^d X > 0, we set 

D\ ■= {xo,yQ) + XDi 

and 

dDy^{x,y) = d\st[{x,y),dDx). 

There exists ci > 0 depending only on Di such that for any (xo,yo) G any A > 0 and all 
f G L°°{Dx), / > 0, the solution z of 

in Dx, 
on dDx 



(3.6) 


satisfies 


z{x,y) 

dD^{x,y) 


> 


Cl A ^ ^ /(x', y')5D^ (x', y') dx'dy', for all (x, y) e Dx- 


Proof of Lemma 5.5, By translation invariance, we may assume xq = yo = 0. If z solves (3.6) in 


Dx, then Z{X,Y) := z{XX,XY) solves 

{-AZ = fx{X,Y)-.= X‘^f{\X,XY) in Hi, 

I Z = 0 on dDi. 
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The inequality for A = 1 is well known; see [9]. Using the fact that 

( 5 di (A"^x, X~^y) = (x, y) for all (x, y) £ Dx, 

and changing variables, it follows that 


(3.7) 


z{x,y) = Z{\-^x,\-^y) > ci5dA>^~ x,\-^y) JJ^ /a(X', U') dX'dU' 

= ci6D,{x,y) [[ f{\X\XY')5D,{XX\\Y')dX'dY' 


Di 


= ciA 6 D^{x,y) f{x ,y)5D;,{x ,y)dx dy , 


Dx 


which proves the lemma. 


□ 


Proof of Lemma 3.2. Starting from the lower estimate (2.14) on {x = 0}, i.e. 

w(0, y) > Cp y^~^ -Cy for 0 < y < p, (3.8) 

the proof is done in three steps (cf. Fig.2 below). 

Step 1. Lower estimate of v in the region Wp n {y > We claim that there exist 

constants K > 0 and p £ (0, L) (depending on v) such that 


j{x,y)>^y^ ^ for (x, y) £ Wp n {y > 


(3.9) 


Let p be given by Proposition 2.3 Using the lower estimate (3.8) on {x = 0}, the fact that 


i;3;(0,y) = 0 and Vxx > —C (cf. (2.18) and (2.4)) and Taylor’s expansion, we obtain 

Op. 


hence 


v{x, y) >Cpy^ - Cy - C x'^ for (x, y) £ uop 
(x, y) > (cp - CK^~^ - Cy^) y^~^ for (x, y) £ Wp n {y > Kx‘^/^^~^'>}. 


The claim (3.9) follows by taking p < {cpfACY^^ and K > {cpjAC) 

Step 2. Harnack-type estimate in suitable boxes near the boundary. We claim that there exist 
constants c, c > 0 and p £ (0, L) such that, for all x £ (0, p/2) and all A £ (0, x/4). 


v{x,y) 


> cA^-^P 


' x—X J 0 


\vy \ dx dy ) — cA for 0 < 2 / < A. 


(3.10) 


By (2.4) and (2.6), reducing p if necessary, we may assume that 

—Ax > |Vx|^ — C > 0 in oOp. 


Let Di be a domain such that 

(-1,1) X (0,2) C Ill C (-2,2) X (0,2). 

For given x £ (0, p/2) and A £ (0, x/4), we set 

D = Dx^\ := (x, 0) + XDi C (x/2, 3x/2) x (0, p) C uOp . 

Observe that —Ax > fx,x ■= |Vx|^ — C in D with fx,\ £ L°°{D) and fx,x > 0. Since x > 0, it 
follows from Lemma |3.3| and the maximum principle that, for some constant ci > 0, 

(3.11) 


> ciA jj {\Vv\P{x',y') - c} dD{x',y')dx'dy'. 


By Holder’s inequality, we have 


D 


\vy\ dx dy < 


Ivyl^doix ,y)dx dy' 


D 


D 


^{x\y)dx'dy 


P-1 

P 
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Using (3.7), we see that 

[[ y) dx'dy' = X^ [[ (AX', AX') dX'dY', = [[ 6~^ (X', Y') dX'dY', 

JJd JJdi JJdi ^ 

where the integral on the RHS is finite due to l/(p — 1) < 1 (see e.g. [33|). Therefore, 

JJ \vy\dx'dy'^ < C JJ \Xv\^6d{x', y') dx'dy'. 

Using also Sd^x', y') dx'dy' = CA^, we deduce from (^3.11^ that 


vix,y) 

dD{x,y) 


> cX^-^P 


\vy \ dx'dy' ) — c A . 


D 


Since 5D{x,y) = y for 0 < y < A, the claim (3.10) follows. 

Step 3. Conclusion. Fix x G (0,p/2) (the case x G (—/o/2,0) can be treated similarly). We 
proceed to estimate from below the integral in (3.10). To this end, we choose 

A = 


where K is from (3.9). Note that this implies A G (0,x/4), taking a smaller p if necessary. By (3.9), 
we have 

rx p 2X px p2\ px „ 

/ / \vy\dx'dy'> / / Vydx'dy' = / v{x',2X)dx' > —X{2Xy~^ = cX‘^~^. 

Jx—xJo Jx—xJo J x—X ^ 

Combining this with (3.10) and using pj3 = /3 + 1 we obtain 

v{x,y) 


y 


>cX ^ — cX for0<y<A, 


which implies, reducing p > 0 again if necessary, 

forO<y<A. 

y 2 

Since 2/3/(l — j3) = 2/{p — 2), letting y —)• 0 we get (3.5). 


(3.12) 

□ 



Fig. 2: The scheme of the proof of Lemma 3.2. 


The marks A, B, C eorrespond to estimates (3.8), (3.9) and (3.12), respectively 


The proof of Theorem 3.1 is then completed by the following lemma. 
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Lemma 3.4. Under the assumptions of Theorem 3.1, for each e G (0,1), there exists constants 
C > 0 and p G (0, L/2) such that 

Vy{x,y)>[{vyf~'P{{l + e)x,U) + {p-l)y]~'^-C mwp\{(0,0)}. (3.13) 

Proof. Step 1. Let po be given by ( |2.7| ) and take p = min(/9o, 1) < Ll2. We first claim that there 
exists a constant ^ > 0 such that 

I-/3 


Vy{x,y)>[{vyYP{x,t)) + {p-l)y] v^xdy-A in Wp \ {(0,0)}. 

Jo 


From (2.2), (2.5), (2.4) and ( |2.7| ), we know that 

\ut{-,T)\<C, \vx\<C, Uy > 0, Vxx>-C in uip \ {(0,0)} 

for some constant C > 1. Let A = 3C and set 

ry 

z{x,y)-.= Vy+ / Vxxdy + A. 

Jo 

Observe that, for all (a:, y) Gujp\ {(0,0)} 

-Zy = -Av < \Vv\P + C. 

Since y < p <1, using (3.15), we see that 


(3.14) 

(3.15) 


ry 

z = Vy + I Vxx dy^ + 2C 'A Vy -\- 2C. 


(3.16) 

(3.17) 


Therefore, by (3.16 )-(3.17), we obtain 


< [{Vyf + + C<[vy + C]P + CP< [Vy + 2C]P < zP, 


hence 


{z^ P)y = -{p - l)zyZ P <p-l in (Vp \ {(0, 0)}. 

After integration, it follows that 

z^~P{x, y) < z^~P{x, 0) + (p - l)y < {vyYP{x, 0) + (p - l)y in uJp \ {(0, 0)}, 
hence the claim. 


ex 


Step 2. We may assume 0 < x < p/2 without loss of generality. To prove ( |3.13[ ), we now take 

we get 

1 ry 

[(vy) ^(• 5 , 0 ) + (p-l)y] ^ds- — ' 

L, so t 
y I r(l+£)* 


the average of inequality (3.14) and we get 

l'(l+s)x Ul+e)x 

Vy{s,y)ds > — 


ex 


Vxx dsdy — A . 


By (2.21), the function x i— Vy{x,0) is nonincreasing for 0 < x < L, so that we obtain 

I r{i+£)x 

Jx ^ Jo 

\-h r<y 


j{s,y)ds > [{vyY P{{1 + e)x,0) + {p - l)y] Vxxdsdy-A 


(3.18) 


> [(^'s/)^ ^((1 + 0) + (p - l)y] ^ -C^- A, 


where we used (2.5) in the last inequality. 

Now we recall from Proposition 2.5 that, for some constant Ci > 0 we have 

(vy - Vx - Cix)x <0 mY. 

Therefore, for some constant C 2 > 0, we get 

2 I‘ii+e)x 

Vy{x,y)> — / Vy{s,y)ds -/ {vxis,y) + Cis)ds + Vxix,y) + Cix 

J X ^X J X 


> 


I rU+e)x 

/ Vy{s,y)ds - C2X , 

J X 


ex 
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where again we also used ( |2.5| ). Jointly with inequality (3.18), we conclude that 

Vy{x, y) > [(uy)^“P((l + £)x, 0) + (p - l)y] - Ce 

i.e., &J3]. 


□ 


Finally, using Lemma 3.4 with e = 1, combined with Lemma 3.2, immediately yields (3.2). The 
proof of Theorem 3.1 is concluded. 


4. Regularizing barriers 


The following lemma shows that a suitable local Holder bound of exponent 1 —/3, near a boundary 
point, actually guarantees a bound for the normal derivative at this point. 

Lemma 4.1. Let p > 2, r,d ^ (0,1), d < L, to G [0,r) and xq be such that [xq — r, xq + r] C 
[-L/2,L/2]. Let 

D = {xq- r, xo + r) X (0, d). 

There exist constants Cq = Co{p) > 0 and rjo = Vo{p,T) £ (0; 1) with the follo wing property. Let 
p G ( 0 , 770 ) and u G x (0,T)) be a nonnegative classical solution of (1.1) in Qt, with u = 0 

on T-r. If 


u{x,y,t) < Cpy^ ^ in Dx[to,T), 


(4.1) 


with Cp = (1 — /3) ^dp, K = Cop^ ^(r^ + T — to) and 
u{x,d,t)<Cp {d +p (t — to)^ r^) 

then 


2 1 

-p^-h {t-to)r‘^ -Ay in [xo-r,xo + r] X [to,T), 


i{xo,0,t) < dpp ^{{t-to)r'^) forte(to,T)- 


(4.2) 

(4.3) 

Remark 4.1. Lemma [4.1| is an improvement with respect to m Lemma 2.2], where a similar result 
was proved, but required a small constant instead of Cp in assumption ( |4.1[ ). This improvement is 
crucial in order to obtain the exact power in the upper estimate of the GBU profile in the next 
section. We note that assumption (4.1) is essentially sharp, since u = Cpy^~^ is an exact solution 
of (1.1) with u = 0 and Uy = 00 at y = 0. 

Proof. Let us define the comparison function 

y2 _ 

z = z{x, y, t) = Cp [{y + (p{x, t))^~^ - p^~^ix, t)] - in T» x [to, T), 


(4.4) 


with 




where 77 > 0 and a G (0,1). Let us denote by C possibly different constants only depending on p 
(often through the value of /3). We first notice that there exists C > 0 such that 

0 < (fit < Cp^~i^r^ipi^ , \Tx\‘^ < Cp^~^{t-to)(p^"^f^ , \Txx\ <Cp^~^{t-to)p'^ . ( 4 . 5 ) 

Moreover, if a is sufficiently small (depending only on p,T), we have 

y{y + (/7)^A < a (1 + ) < Cp(l - jl) = dp . ( 4 . 6 ) 

This implies 

i^yjy + '' ^ 

d p j 


\zy?‘ = dl{y + p) (^l 


) <dl{y + p) 
V 
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and since 

we deduce 


\zx\‘^ = dl (^{y + <f) 

\Vz\P = {zl + zl)^ < (E {y + ip)-P^ [l + (y + 


Since, from (4.5), we have < Crj^ — tQ)p^ by taking rjo = r]Q{p,T) sufficiently small, 

it follows that {y + < (j_ Therefore we have 

\Vz\P <dP{y + p)-Pd l + dy + pf^plp-"^^ , 

which implies 

iV^r < (y + p)-P^ + C[y + pf-^plp-^d 
<dP{y + p)-P^ + Cplp-d-^. 

Thus, for (x,t) G D X {tQ,T), we estimate 

zt-Az-\Vz\P > dp[{y + p)~^ - p~d]{pt - 

+dp/3[iy + p)-d-^-p-f^-^]pl 

+dpP{y + p)~^~^ + K 

-dp (y + p)-P^ - Cplp-d-^. 

Using 

dp/3(y + p)~^~^ = dFp{y + p)~P^ 

we deduce 


zt 


- Az - \Vz\P > -Cp ^{\pt\ + \pxx\)) - Cp d VI + 


and thanks to (4.5) we conclude 

zt-Az- iV.zl?’ > K - CV + T- to)V"^ ■ 

In particular, we have 

zt-Az> \Vz\P inDx (to, T) (4.7) 

with the choice k = C{r^ + T — to)y^“^ (which in turn guarantees (4.6) for rjQ = r](}{p,T) small). 

On D X {to}, as well as on the lateral boundary part {xq — r, xo + t} x [0, d] x (to, T), we have 
p = 0 and so 

7/2 

1-/3 y \ 

z = Cpy P - K— > u 


by (4.1). Next, on the part [xq — r,xo + r] x {0} x (to, T), we have, for to < t < T, 

-E) = zi-,-E) = 0. 


On the remaining part [xo — r, xo + r] x (dj x (to,T), thanks to (4.2) and to the fact that the 
expression in (4.4) is a decreasing function of p, we have 

r 12 1 d^ 

u{x,d,t)<Cp {d + p{t — to)^ — ri^~d {t — to)r^ — k—< z{x,d,t) ■ 

By the comparison principle, we deduce 

u< z in T) X [to,T). (4.8) 

In particular, we have 

Uy(xo,0,t) < Zy{xo,0,t) 

hence (4.3). □ 
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5. Proof of main result: the upper estimate 


We first establish the following 


Proposition 5.1. Let p G (2,3], q > p — 1 and let u be as in Theorem 1.1. There exist k G (0,1) 
and xi,yi,a > 0 such that we have 

+ + y)u'^ <0 in (0, xi) x (0, yi) x {T — a, T). 


Ucr 


Remark 5.1. It follows from Proposition |5.1| and estimate (2.14) that there exists a constant c > 0 
such that, for y > 0 small. 


Uxx{0,y,T) = lim 
x^O 


Uxix,y,T) 


<-fc(l + y)(^ ^1-/3 ) 


—c < 0. 


Proof of Proposition \57l\ It is divided into several steps. 

Step 1. Preparations. We consider the auxiliary function 


J(x, y, t) = Ux + c{x)d{y)F{u) in D x (0, T), 

where D = (0, xi) x (0, yi) and the smooth positive functions c, d, F will be chosen below. Our aim 
is to use the maximum principle to prove that, for sufficiently small xi, yi, <7 > 0, there holds 


J<0 in H X [T — cr, T). (5-1) 

This will be done in the subsequent steps. 

Step 2. Derivation of a parabolic inequality for J. The following basic computation was 
made in m For completeness and for the convenience of readers, we repeat it here. We have 

Jt — Uxt T cdF Ut, 

AJ = Aux + cdF'Au + cdF'^lVul"^ + 2c'dF'ux + 2cd'F'uy + {cd" + c''d)F. 

Then we obtain 


Jt-AJ = {\Vu\P)x + cdF'\Vu\P - 2c'dF'ux - 2cd'F'uy - cdF"\Vu\^ - {cd" + c"d)F. 
Using Ux = J — cdF, we write 

(|Vu|P),, = p\Vu\P-^Vu-Vux 

= p\Vu\P-‘^Vu • VJ - p\Vu\P-^Vu • V{cdF) 

= p\Vu\P-‘^Vu • VJ - p\Vu\P-‘^{cdF'\Vu\‘^ + Uxc'dF + Uycd'F), 

hence 



(iVnD. = 

p\Vu\P-^Vu ■ VJ - pcdF'\Vu\P 

-pc'dF\Vu\P-^J + pccd^F‘^\Vu\P-‘^ -pcd'F\Vu\P-‘^Uy. 


We also have 


-2c'dF'ux = -2c'dF'J + 2cc'd'^FF'. 


So we get 

II 

< 

1 

aJ + b-VJ 

-{p- l)cdF'\Vu\P + pcc'd‘^F^\Vu\P-^ - pcd'F\Vu\P-^Uy 
+2c(fd?FF' - 2cd!F'uy - cdF''\Vu\^ - {cd” + c"d)F, 

(5.2) 

where 

a = 

-p(fdF\Vu\P-^ - 2c dF' and b = p\Vu\P-‘^Vu. 

(5.3) 

Let 


VJ = Jt- AJ -aJ -b- VJ. 
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We can rewrite the above equality as follows 

VJ 


cdF 


F' F" r” d" 

= -(p-i)-ivur--iv.p---- 

—2 — —Uy — p—\Vu\^~‘^Uy + pcdF\Vu\^~‘^ + 2c'dF' 
d F d 


(5.4) 


Step 3. Estimation of the RHS of (5.4). We now specialize the previous computation to 
the following choice: 


F{u) = 

d{y) = y~'^p{y), 

c{x) = kx, 

where /c G (0,1), (; > 1, 7 > 0 and (/? > 0 is a smooth function with ip' > 0. Using 

d'{y) = -72/“’^" V(2/) + y~'^‘f'{y), 
d”{y) = 7(7 + - 27y"'^“V'(2/) + y~^y^"{y), 

the equality (5.4) implies 

y^VJ 


(5.5) 

(5.6) 

(5.7) 


cdF 


y‘^\Vu\p ! i/2|v^|2 ^ yy^' 2/V" 

= -(P-I)?-9(9-1)-2-7(7 + 1)+ 27- 

U U‘^ ip p 


(5.8) 


+2 q(^+ p(^—^^'^y\VuY’ “^Uy+pky^ '^pu'^\Vu\^ + 2kqv!^ ^y^ '^p. 

Also, taking u,xi,yi G (0,1) sufficiently small and setting Q = (0,xi) x (0,yi) x (T — cr,T), we 
have, by ( 1 ^, (( 2 ;§, 

\ux\ < Cx in Q, 

as well as 

y <u < Cpy^~^ + Cy'^, I < Uy < dpy~^ + Cy in Q. (5.9) 

In particular we have, close enough to the singularity. 


y\Vu\P ^Uy = y{uyY ^ 


1 + — 


u. 


( p - 2)/2 


< y{ui 


\p-i 


1 + 


p-2 


2 Vu 


= y{uyY ^ + ^^-i^y{uyY ^{uxf 
^yiuyY-^ + Cy'^x^ 

with m = min(l, 2/3). In particular, for p < 3 (i.e. /3 > |), we have 

y\Vu\P~‘^Uy < y{uyY~^ + Cyx‘^. 
Similarly, using (5.9) we estimate 


+ 2kqu'i-^y‘^-^p < kqCpy^'^-^^^^-d^+^-'^ < qCpy^+^+'^^^-f^^-'^ 
et 

y^iUyY 


for any k G (0,1). Consequently, we get from (5.8) 

y‘^VJ 


cdF 


< -{p-l)q- 


u 


7(7 + 1) 


+ 29(7 


u 


P-1 


p J u 

+ 27 ^ _ ^ + Cy{x‘^ + qpy^+'i^^-^^-Y 

p p 

We will conclude Step 3 through the following lemma. 


(5.10) 


in Q. 
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Lemma 5.2. Let p G (2,3], g > p — 1, and take 7 = g(l — (5) and ip{y) = 1 + y in (5.6). There 
exist xi,yi, cr > 0 sufficiently small such that, for any k G (0,1), we have 

VJ <0 in Q = (0, xi) X (0, yi) x {T — a,T). 


Proof of Lemma 5.2. To shorten notations, we set 

^ y,t) = ^ >0, 0 = 0 (x, y, t) = y{uy)P~^ > 0 , V’ = V’(y) = 

Notice that '4>{y) = is small provided yi is sufficiently small. 

We wish to show that the right-hand side of (5.10) is nonpositive in Q. To this purpose we 
distinguish two cases according to whether ^ < 1 — /3 or not. 

Case 1: f, < 1 — ffi 
By Young’s inequality, we have 

y{uy)P 


p 7 


u, 




+ 7 - 


P u 


p -1 




qP-'-yP 


(5.11) 


Using (5.9), and 1 — /3 = (p — 2)/(p — 1), we have 


p uP ^ ^ p/ 'ijj{y)Y{cpy^ f^ + Cy‘^)P ^ 

-li,p-2 — "1^ I 


qP-ffiP- 


7 


< 




qP-^yP-'^ 

\l + Cy^^y 


p-i 


The precise value of Cp in (2.11) and the choice 7 = g(l — /?) then imply 


7 


gP^^ ^ 7 / 3 ( 1 + ' <g(l-/?)/?-p/3^(y) + Cyi+^ 


for some C (possibly depending on q). Hence from (5.11) we obtain 


p 7 


U; 


,)P-^ < (p - l)qy^^ + q{l - /3)(3 - p^y) + Cy^^^ . 


Therefore, using ^p" = 0, (5.10) implies 

y‘^VJ 


cdF 


< -g(9-l)e' + 2(?(7- 




- 7(7 + 1 ) + 9(1 - /3)/3 - P/3V’(y) + 27 V’(y) + Cy{x^ + y^). 

Now we remark that the function 


-q{q-l)e + ‘^q(l- 

is increasing for ^ < 1 — ,0 and y sufficiently small, so we get 




y^VJ 


< -q{q - 1)(1 - / 3 )" + 2 ( 7(7 - ^{V)) (1 - / 3 ) 

- 7(7 + 1) + 9(1 - / 3)/0 - + 27V’(y) + Cy { x ^ + y ^) 

which implies, using 7 = (1 — /3)g, 

< -Pl^ifiy) + Cy{x‘^ + y^) = -p/3j^ + Cy{x^ + y^). 
cdr 1 + y 

Therefore, we have VJ < 0 provided xi, yi are taken sufficiently small. 

Case 2: | > 1 — /3. 
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With the above notations, and using ip" = 0, (5.10) can be written as 

y^VJ 


cdF 


hence 


where 


<-{p-l)q^Q-q{q-l)^ - 7(7 + 1) 

+ 2^(7 - V’(y))? + p(7 - ^( 2 /))© 

+ + Cy{x^ + /) 

ii^V T 

< i?(y,^,0) + C'y(x^ + /), 


(5.12) 


H{y,^,Q) = {p[qil -/3) - i^iy)] -g(p-l)^}0 

- q{q - 1)^2 + 2 q[q{l -(3)- ^|;iy)]C (5-13) 

- q{l - /3){q{l -/?) + !)+ 2g(l - /3)iP{y ). 

Subcase 2.1: p[q{l — /3) -i’iy)] -9(p-1)? < 0. 

Using ^ > 0 and 0 < ip{y) < yi in Q, we first observe that this implies H < q f{^), where 
/(O := -{q - 1)?' + 2g(l - /3)^ - (1 - /3){q{l -/?) + !)+ 2(1 - f3)y,. 

Computing the reduced discriminant of this trinomial we notice that 

A = q\l - Pf -{q- 1)(1 - I3){q{l - /3) + l) + 2(g - 1)(1 - (3)y^ 

= {1 - /3)[l - q/3 + 2{q - l)?/i] < 0 

provided q > ^ = p—1 and yi is sufficiently small. Therefore /(^) < 0 for every ^ > 0 and we have 


p[q{l - (5)- V’(y)] - q{p - 1 )^ < 0 


H < q max/(,^) < 0. 
~ ?>o ^ ^ 


Hence VJ < 0 provided yi is sufficiently small. 

Subcase 2.2: p[q{l — l3)-'il^{y)] -g(p- 1 )? > 0 . 

Here, we use ( |5.9[ ) to estimate 

0 < [cp(l - /3W-\l + = /3(1 + < /3 + Cy^+^ 

Since p[q{l — (3) — "0(2/)] ~ q{p — 1)^ > 0; h follows that 

H < {p[q{l -13)- i’iy)] - q{p - 1)^} 13 
- q{q - 1 )?^ + 2g(g(i -13)- V’(y))^ 

-q{l-13) (g(l -f3) + l)+ 2q{l - /3)i;{y) + Cy^+f^ 
which yields, using pf3 = 1 + f3, 

H<q{e-C + l3il-l3)-q[^-{l-f3)]‘'-2i;{y)[C-il-f3)]] 
-(l + f3)i;iy) + Cy^+f^. 

Now we observe that, for any y G [0,yi], the function 

</>(0 := e' - e + /3(1 - /3) - q[i - (1 - /?)]' - 2i^{y)[i - (1 - /?)] 


is concave and we have 


</)'(!-/3) = 1 - 2/3 - 2iA(y) < 0 


(5.14) 
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since /? > ^ ( p < 3). Therefore, we have 

<^(0 ^ '/’(I “ /5) = 0 for all ^ > 1 — /3 


(5.15) 


and we conclude that 
hence 


H <-{l + (3)i;{y) + Cy^+P, 


,,2'P T 

< -(1 + I3)'il^{y) + Cy^^^ + Cyix"^ + /) 

<-(l + /3)^ + Cy(x2 + /)<0 
1 + y 


□ 


provided xi and yi are taken sufficiently small. The proof of Lemma 5.2 is complete. 

Continuation of proof of Proposition \5.1\ 

Step 4. Initial and boundary conditions for J. First observe that, for each T' < T, we 
have 


u < Cy in D X [0, T'] 
for some C = C{T') > 0. Since j < q, we have in particular 

J eC(Dx [o,T))nc^’\D X (o,r)). 


(5.16) 


(5.17) 


Also in view of (5.16) and 7 < g — 1, the coefficient a{x,y,t) of the operator V (cf. (5.3)) satisfies 

a is bounded in 71 x (0,T') for each T' < T. (5.18) 

Next, since w = Ux satisfies 

wt — Aw = p\S/u\^~^Vu -Vw in Qt, (5.19) 

and is nonnegative nontrivial in ( 0 ,L)^ x ( 0 , T), by the maximum principle and after a time shift, 
we may assume that 


Ux < 0 in (0,L)^ X (0, T). 


(5.20) 


Let now xi,yi,cr be given by Lemma 5.2 and assume a < Tf2 without loss of generality. By 


(2.1), (5.19), (5.20), (1.11) and Hopf’s Lemma, there exist constants ci,C 2 > 0 such that 

Ux<-ciy on {xi} X (0,yi) X (T/4,r), 

Ux<-cix on (0,xi) X {yi} X (r/4,r), 
u<C 2 y on {xi} X (0,yi) X (0,r). 

With the above estimates, we check the function J on the lateral boundary: if y = yi, we have, 
by d^ , 

Jix,yi,t) < [-Cl + kyf'^{I + yi){cpyl~^ + cyi)‘']x <0 on (0,xi) x {yi} x (r/4,r), (5.21) 

if k is sufficient small. If x = xi, we have 

Jixi,y,t) < -ciy + fcxic|y'?"'^(l+ y) 

< [-Cl + A:xiC 2 y?“'^“\l + yi)]y < 0 on {xi} x (0,yi) x (r/4,T), (5.22) 

if k is sufficiently small, since 7 = g(l — (3) < q — 1. Moreover, we clearly have 

J(x,0,t) = 0 on (0,xi) X {0} X (T/4,T), (5.23) 

«7(0,y,t) = 0 on {0} X (0,yi) X (r/4,r). (5.24) 

Finally, we recall that there exists C 3 > 0 such that 

Uxix,y,T - a) <-C 3 xy in D. 


(5.25) 
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(This is a parabolic version of “Serrin’s corner Lemma”; see |27L p.512]). Now (5.25) implies 
J{x,y,T-a) < -csxy + kxcly‘^~^{l+ y) 

< [-C 3 + + yi)]x?/< 0 in (5.26) 

if k is sufficient small, since 7 < g — 1. Then (5.1) follows from Lemma 5.2, (5.21 )-(5.24), (5.26|) 


and the maximum principle. Note that the use of the maximum principle is justified in view of 

□ 


(5.17) and (5.18) (or, alternatively, of the fact that a < 0). 


The proof of Proposition 5.1 is complete. 


By combining Proposition 5.1 and Lemma 4.1, we shall now prove the upper estimate. 

Proof of Theorem 1.1: the upper estimate in (1.8). It suffices to prove it for x > 0 sufficiently small 
(the case x < 0 will follow by considering u{—x,y,t)). By Proposition 5.1, we know now that, for 
some k>0, q>p—1 and to G ( 0 , T), we have 


u, 


+ kxy ^)(l + y)n ^<0 in Q = ( 0 , xi) x ( 0 , yi) x (to, T). 


Integrating over (0, x) and using (2.10), we obtain 

u^"'?(x, y, t)>{q- l)fey + u^"'?( 0 , y, t) 

x 2 


(5.27) 


> (g - + [cpy^-f^ + cyY“'' m Q. 


Starting with this estimate, we shall now apply the regularizing barrier lemma (Lemma 4.1). 


Fix xo G (0,xi). Let rj G (0,7o)) where % is given by Lemma 4.1, and set 

Xq 

r=y, d = r]Xo, D = (xq - r, xq + r) x {0, d). 


(5.28) 


Next we recall k, given by Lemma |4.1[ 

K = Co{p) + T - to). 


(5.29) 


We shall also use the notation t = t — to. 

We claim now that there exists 7 G ( 0 , 70 ) such that, for any xq sufficiently small, we have 


Cpy^ ^ > 0, 


and 


Cp [{d + yr iWr iW )^“^ — >0, 


in D X [to, T) (5.30) 


|cp[(d + yr 1-/3r 1-/3)^ ^ — 7 ^ 


l-q 


in [xo - r,xo + r] X [to,r). (5.31) 


-[cpd^-^ + cd^Y '^ < (g - l)A:^d-9(i-/3) 


Assume for the moment that ( |5.30 )-(5.31 ) hold; together with (5.27), this implies 

„,2 j 1-9 


'‘{x,y,t) > \ Cpy^ ^ 


in D X [to, T) 


and 


“^(x, ct, t) > |cp[(i + r/T 1 -/ 3 r 1-/3)^ ^ — 7^ 


l-q 


in [xo -r,xo + r] x [to,r), 
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and so both (4.1) and (4.2) will hold. We may then apply Lemma 4.1 to deduce 


/{xo,0,t) < dpT] P[Tr‘^) 1-/3 in(to,T). 


(5.32) 


At time t = T this gives 


2/3 


uyixo, 0, T) < xo = Crj-^ x^ """ 


which, jointly with ( 2 . 8 ), implies the upper estimate in ( 1 . 8 ). 


To conclude, we are thus left to prove (5.30)-(5.31). We note right away that the first part of 
( |5.30 ) is true whenever d is small enough, a condition which holds as soon as xq is small enough 
(independently of 77 ). Similarly, observe that if xq - hence d - is sufficiently small (independently 
of r]), then 


7 , 2 , 

1-/3 y 

CpV ^ - Ky 


1-9 


- [cpV^ ^ + cy'^] 


, 211-9 


< {Cpy^-f^r^ [(1 - - (1 + 

in H X {to,T). Here and in the rest of the proof, C denotes a generic constant independent of xq 

and rj. Consequently, (5.30) holds as soon as r] is sufficiently small. 

1 2 


In order to prove (5.31), setting Q = rjT^-i^, we write 

92 


Cp[(d + r/r 1-/3r 1-/3)^ ^ — 77 ^ ^rr^] — Ay = Cpd^ ^ 




> Cpd} ^ 


1+/3 \ 1-/3 


x\ 1-/3 

- —d^+f^ 


2c„ 


/ _j_ K , , 

1 — (r 1-/3r 1-/3 j- (rfr) 

\ / 2c7i 


1+/3 


By ( |5.29[ ), (5.28), it follows that 


Cp[{d + riT^r^)^~^ — r]^~^Tr‘^] — K— > Cpd^~^ 1 — Cxq"*"^ (T + (xg + T)r 7 Q) >0, (5.33) 


provided 


Cxl+^ {T + {xl + T)yl) <1. 


(5.34) 


Note that (5.34) is true for xq sufficiently small, which guarantees the first part of (5.31). 
Next notice that the convexity inequality (a + > a^~'^ + (1 — q)a~^b implies 

[cpd^-^ + cdY“^ > {cpd^~^Y~‘^ - 


Combining this with (5.33) and (5.34), it follows that 


|cp[(d + r7T 1-/3r 1 -/ 3 )^ ^ — 77 ^ ^3‘r^]—Ky| — \cpd^ ^ + cd^] 


211-9 


< {cpd^-^f-^ I [1 - CxJ+^ (T + (xg + T)yl)\^ - 1 [ + Cd-'?(i-^)+2 

< (T + {xl + T)yl) + 

(t + {xl + T)yl + 77'+^) xg d-'?(i-«. 


< 


Therefore, (5.31) is satisfied as soon as rj is chosen sufficiently small. 


□ 
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k:(o) = v"""- 


6 . A HEURISTIC EXPLANATION OF THE SINGULARITY EXPONENTS THROUGH QUASI-STATIONARY 

APPROXIMATION 

A possible heuristic explanation of the appearance of the number 2/(p — 2) in the tangential 
singularity profile ( 1 . 8 ) can be obtained using the idea of quasi-stationary approximation along 
the family of ID steady states. 

Recall the following family of ID steady states, given by the translates of the reference solution 

y{y) = Cpy^~^, 

i.e. 

Va{y) = V{y + a) -V{a), y > 0, o > 0. 

These special solutions verify 

-C = C, K(0) = 0, 

The idea is then to look for an approximate solution obtained by modulating in a, or moving on 
the manifold of steady-states (14)a>o- More precisely, we set U = Uapprox given by 

Uix,y,t) = V{y + h{t,x)) -V{h{t,x)), (6.1) 

which amounts to parametrize the solution by a = h{t,x). In particular, we have Uy{x,y,t) = 
V'{y + h{t,x)). The function h{t,x) is positive for t < T and must satisfy h{T,0) = 0 so that 
Uy {0, 0, T) = oo. 

Note that this Ansatz means in some sense that —Uyy ~ {uy)^ and ut ~ Uxx near the singularity, 
already giving a rough clue to the parabolic nature of the scaling of the profiles in t and x. 

With the above Ansatz, one has an interpretation of the lower estimate of the tangential profile in 
( 1 . 8 ) as being a consequence of the constraint Uxx > ~Ci which comes from the maximum 
principle (cf. Proposition 2 . 1 ). Indeed, Uxx > —C and Ux{0,y,t) = 0 imply that 

Ux > —Cx, X > 0. (6.2) 

For t < T, restricting without loss of generality to x > 0, we note that 

Uxix, y, t) = dp [(y -h h{t, x))~^ - x)] hx > -dphxh~^{t, x), 

where we used hx > 0 due to Ux < 0 , and that 

Ux{x, h{t, x),t) = —chxh~^{t,x). 

Consequently, (6.2) is equivalent to hxh~^{t,x) < Cx. By integration, it follows that h^~^{t,x) < 
0) -|- Cx^. Letting t we get h{T, x) < Cx^l^^~^\ which leads to 

Uy{x, 0, T) = V'(h{T, x)) = dp(h{T, x))"^ > = Cx-2/(p-2). (6.3) 

The fact that the upper estimate in (1.8) is exactly of this order means that the constraint Uxx > — C* 
is satisfied in a minimal way by the parabolic flow. 

The same analysis can be done with the time rate as well and actually enables one to recover 
also the exponent l/(p — 2 ) of the time rat^ namely 

||Vu(t)||oo~(T-t)-'/(^-") 

(the lower estimate is always true - see |11( IT 8 l l32] - whereas the upper estimate is only known for 
monotone increasing solutions in ID; see |18] and cf. also [32|). This time the essential constraint is 
\Ut\<C (cf. Proposition |2.1| ). Indeed, one can easily see that \Ut\ < C is equivalent to \ {V{h))t\ < 
C. Since h{T, 0) = 0, we thus have V{h{t, 0)) < C{T — t), i.e. h{t, 0) < C{T — or 

Uy{0, 0, t) = V'{h{t, 0)) = dp{h{t, 0))-^ > C{T - t)-/3/(i-/3) = c{T - t)-i/(P-2). (6.4) 


■^this observation doesn’t seem to have been made in previous work on the time rate. 
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Since the rates (6.3) and (6.4) violate the self-similar structure, or natural scaling, of the equation 


(cf. Remark (a) in Section 1.3), so one can say that the maximum principle here wins against 
self-similarity. 


Remark 6.1. For p > 3, the proof of the upper estimate in (1.8) fails at the level of inequalities 
(5.14)-(5.15). Actually, it can be seen along the lines of the proof of Lemma 5.2 that VJ > 0 
in some regions near the singularity (more precisely, where yUyU~^ ~ (1 — /5)_|_ and yuy~ ~ /3). 
However this might be technical, and it is presently open whether or not the actual behavior of u 
changes for p > 3. 

As for the above heuristic argument, although it does not a priori make a difference between 
the ranges 2 < p < 3 or p > 3, it is not clear if such an argument can suggest more than a lower 
estimate of the profile. Indeed, we stress that the heuristic argument gives a justification of the 


lower estimate in (1.8) only in view of the one-sided estimate 

u-TT > —C. 


(6.5) 

For p < 3, our proof of the upper estimate also shows that the estimate Uxx > — C* is really optimal 
and that Uxx is discontinuous near the singularity at t = T (cf. Remark 5.1). If one could show for 
p > 3 that Uxx remains continuous (i.e., has a zero limit) near the singularity at t = T, then the 
proof of Theorem |3.1| (as well as the heuristic argument), would imply that the final profile of Uy 


is more singular than (1.8). 


However, such property of Uxx should be rather unstable and proving this might be quite delicate. 
Indeed, for any p > 2 and any a > (p — l)/(p — 2 ), a simple computation shows that the function 

u{x, y, t) = Cp [(|x|2“ + iT- t)" -h y)^~^ - (|xp“ + {T - t)“)^"^], 

modeled after ( |6.1[ ), is a solution of 

ut — Au = \Vu\^ + f 

in Q = (—1,1) X (0, 1) X (0, T) with u{x, 0, t) = 0 and some / G L°°{Q). Moreover u has an isolated 
gradient blowup point at (0, 0 , T), u satisfies (6.5) and \ut\ < C, but Uxx is continuous near (0,0,T) 
if and only if a > (p — l)/(p — 2). 
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